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ABSTRACT 

The formation of thick stellar disks in spiral galaxies is studied. Simulations of gas-rich young galax- 
ies show formation of internal clumps by gravitational instabilities, clump coalescence into a bulge, 
and disk thickening by strong stellar scattering. The bulge and thick disks of modern galaxies may 
form this way. Simulations of minor mergers make thick disks too, but there is an important difference. 
Thick disks made by internal processes have a constant scale height with galactocentric radius, but 
thick disks made by mergers flare. The difference arises because in the first case, perpendicular forcing 
and disk-gravity resistance are both proportional to the disk column density, so the resulting scale 
height is independent of this density. In the case of mergers, perpendicular forcing is independent of 
the column density and the low density regions get thicker; the resulting flaring is inconsistent with 
observations. Late-stage gas accretion and thin disk growth are shown to preserve the constant scale 
heights of thick disks formed by internal evolution. These results reinforce the idea that disk galaxies 
accrete most of their mass smoothly and acquire their structure by internal processes, in particular 
through turbulent and clumpy phases at high redshift. 

Subject headings: ISM: structure — Galaxy: disk — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

In the cosmological model of galaxy formation, gas-rich 
mergers are necessary to preserve d isks and build spiral 
galaxies (e.g.. lRobertson et al.ll2006ft . These mergers also 
scatter disk stars and produce the thick disk component 
(jQuinn et all 119931 : IWalker et ai]|1996fh There is how- 
ever increasing evidence that most of the baryonic mass 
enters in cold flows, which are smoother and less dis- 
ruptive than mergers (lAgertz et all 120091 : lOcvirk et al.1 
l2008tlKeres et al.M2009h . Cold flow accretion is ex pected 
to be most active at z ~ 2 (|Dekel et all 12009ft . Ob- 
servations of galaxies around that epoch show a large 
popul ation of rotating disks that are clumpy and turbu- 
lent (lElmegreen fc Elmegreer][2fM iGenzel et all [2008: 
IShapiro et alj 120081 : lElmegreen et all l2009allbl ). The 
clumps in these galaxies have masses up to a few 10 s M , 
sizes of around a kpc, an d likely formed intern ally by 
gravitational instabilities (|Elmegreen et^l200l . Cold 
flows have an advantage over mergers as an explanation 
for high-redshift clumpy galaxies (Agertz et al. 2009: 
Ceverino, Dekel & Bournaud 2009). This is because the 
formation of the observed clumps requires a relatively 
high disk velocity dispersion for the ambient Jeans mass 
to be large. Then a high fraction of the baryons has 
to be in the primordial disk to keep the system unsta- 
ble: only relatively smooth gas accretion can do that 
(|Bournaud fc Elmegreer]|2009t ' iDekel et al.ll2009t) . 

The origin of thick disks then comes into question. If 
most of the mass assembly is from smooth gas flows and 
rapid internal evolution, then thick disks have to form by 



internal processes, such as stellar scattering off of clumps 
and disk hea ting by gravitationa l instabilities that 
make clumps (|Bournaud et all 120071 hereafter BEE07). 
Clum py disk evoluti on can form classical bulges (iNoguchil 
119991 : llmmeli et all 12004 lElmegreen eTalTl2008f ). and 
exponential disks (Bournaud, Elmegreen & Elmegreen 
2007, hereafter BEE07), but there are no distinct fea- 
tures in these remnants that exclusively point to clump- 
driven evolution. Perhaps thick disk formation by clump 
scattering is a more revealing remnant in modern galax- 
ies. 

Here we explore the properties of thick disks formed 
internally in unstable, gas-rich, clumpy disks. We com- 
pare with simulations of merger-induced disk thickening. 
We find that thick disks formed internally have a nearly 
constant scale height with radius, as observed for real 
galaxies. This is unlike the case for merger-induced thick 
disks, which always flare. We also investigate how much 
the thick disk changes during later gas accretion when 
the thin disk forms, and find that the scale height de- 
creases, but remains constant with radius. 

2. PROPERTIES OF OBSERVED THICK DISKS 

The primary remnants of the earliest stage of galaxy 
formation are the bulge and thick disk components. The 
thick disk in the Milky Way conta i ns no stars younge r 
than 8 Gyr (|Gilmore et "all 119851: Heddv et all l2006h. 
Thick disks are u biquitous (jDalcanton fc Bernsteinll2002l : 
ISeth et all 120051 ) . and contain 10% to 25% of the bary- 
onic mass of spiral galaxies. In the disk alone, thick- 
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Fig. 1. — Edge-on view of the Milky Way analogue NGC 891 
in the i?-Band from 2Mass data (a). Vertical luminosity cuts in 
the 3 < r < 4 kpc (b) and 6 < r < 8 kpc (c) radial bins, avoid- 
ing any foreground star, show the thin mid-plane dust lane, the 
thin stellar disk, and the thick stellar disk. The thick disk has a 
relatively constant scale-height. Even the deepest HST/ACS data 
for NGC 891 can be fitted with a thin +thick disk model having 
uniform scale-height (Ibata ct al. 2009). 

to-thin mass ratios range from 10% to 50% (Yoachim 
& Dalcanton 2 006, hereafter YD06); i n the Milky Way 
it is 20 - 25% (IGilmore fe ReidlH983l : IChen et al.l[200ll 
iRobin et al.ll2003f) . 

The structural properties of thick disks that are rele- 
vant for comparisons with our models are: 

• The thick-disk sech 2 scale height, zo, is gener- 
ally between 1.4 and 2.8 kpc, with a large scat- 
ter (Dalcanton & Bernstein 2002, YD06). Equiva- 
lent exponential-fit scale heights, h z , are about half 
these values. 

• Zo is about constant with rad ius (YD06). In the 
Milky Way (Robi net al.l I2003T ) and other galaxies 
fe.g. Jlbata et al.l | 2009 !), h z also does not have sig- 
nificant variations until R25 and beyond. This con- 
stancy of h z is illustrated here in Figure [TJ where 
vertical profiles from the inner and outer parts of 
NGC 891 are presented. A constant h z implies that 
thick disks viewed edge-on have disky isophotes. 
Signs of boxiness were sought but not found (Dal- 
canton & Bernstein 2002). 

• Thick and thin disks have relatively similar radial 
scalelengths for most galaxies (YD06). 

Bulges and thick disks have enhanced a elements com- 
pared to Fe (e.g., iLecureur et alj l2007t iZoccali et al.l 
120071 ). suggesting a star formation timescale shorter than 
1 Gyr. Other chemical similarities between the bulge 
and the thick d i sk of the Milky Way were discussed by 
iMelendez etHI (|2008l ). 

3. SIMULATIONS 

3.1. Unstable disks at high-redshift and thin disk growth 

The simulations used here are similar to those in 
BEE07, but with an increased resolution to see the verti- 
cal disk structure. The iV-body particle-mesh code uses 
a sticky particle scheme to model cold gas dynamics. The 
spatial resolution (softening size) is 35 pc. Two million 
particles are used for each component of gas, stars, and 
dark matter. We ran 3 models: 



• Model 1 starts with a stellar mass of 2xl0 10 M Q 
and a gas fraction of 60%. 

• Model 2 starts with a stellar mass of 3xl0 10 M Q 
and a gas fraction of 50%. 

• Model 3 starts with a stellar mass of 4xl0 10 M© 
and a gas fraction of 40%. 

The models start with idealized, uniform disks having 
masses that are plausible for z~2 progenitors of Milky 
Way- like spirals. They rapidly acquire a phase space 
distribution (morphology and kinematics) that is similar 
to that in cosmological simulations of clumpy disk for- 
mation (Agertz et al. 2009; Ceverino et al. 2009) and 
also consistent with that of observed clumpy galaxies (see 
BEE07 and Bournaud et al. 2008). Thus the present 
models should be adequate for studies of disk evolution 
during and after the clumpy phase. 

The initial disks have constant surface density for both 
gas and stars. The external radius is 6 kpc and the initial 
scale height is h z = 500 pc. The ratio of dark matter to 
baryon mass is 1:2 inside the initial disk radius. The 
models are first run for 1 Gyr in full isolation in order to 
identify the effects of internal evolution. 

After this time, additional gas is allowed to accrete 
for another 6 Gyr. This is supposed to model the 
later growth of the thin disk. The gas accretion rate 
is 10 M Q yr -1 , which increases the initial galaxy mass 
by a factor 2 to 3. Accreted particles are positioned at 
a large radius (15 kpc), close to the disk plane: they are 
initially given a vertical distribution with an exponential 
profile and a scale height of 1.5 kpc, and later-on settle 
in a thinner disk. Accreted particles are not all present 
in the simulation initially, they are continuously added 
at the r=15 kpc boundary to model continuous accretio n 
down to low redshifts, as in lBournaud fe Combed (2002). 
Each is given a rotational velocity such that it has the 
same angular momentum as an average particle in the 
thick disk at 1 Gyr. This similarity of angular momenta 
gives the thick and thin disks similar radial scale lengths. 
Thus the additional gas is accreted from a thick layer in 
these models, but is relatively stable and naturally settles 
to a thin disk. 

3.2. Minor mergers 

Many simulations of min o r mergers have been pub- 
lished dWalker et al.l \l99(k iVillalobos fe Helmil 120081 : 
iRead et al.ll2009h . in addition to cosmological simula- 
tion s with with numer ous satellites accreting and merg- 
ing (|Abadi et al.ll2003t ). We performed new merger mod- 
els to compare with the internal evolution using the same 
code and assumptions. We ran four idealized simulations 
of single minor mergers and a cosmological simulation 
where several minor mergers occur with mass ratios and 
orbits prescribed by a A-CDM model. 

The four individual minor mergers were chosen to have 
a mass ratio of 10:1. The primary galaxy is modeled with 
one million particles for each component (gas, stars, dark 
matter). Its stellar mass is 4xl0 10 M and its gas mass 
is 1 x 10 10 M . It starts with an exponential disk of radial 
scale-length 4 kpc, truncated at 15 kpc, and an exponen- 
tial scale-height of 500 pc. 20% of the stellar mass is in 
a central bulge with a Hernquist profile of scale-length 
500 pc. The dark matter halo has a Hernquist profile 
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Fig. 2. — Model 2 of a gas-rich gravitationally-unstable disk gen- 
erating giant clumps of gas and stars, similar to clumpy galaxies 
at z > 1 (see also BEE07). These edge-on views of the gas+star 
mass density show rapid thickening of the disk by the clumps. The 
clumps finally coalesce into a central bulge within a Gyr. Subse- 
quent gas infall forms a stable thin disk inside the thick disk. The 
thick disk scale-height decreases a little when this happens, but the 
separate thick disk component remains. 

with a scale-length of 7 kpc. The dark matter-to-baryon 
mass ratio is 1:2 inside the disk radius. The companion 
galaxy has all masses divided by a factor of 10, and all 
sizes divided by 3.5. The orbital inclination with respect 
to the target disk plane is 35 degrees. The impact pa- 
rameter is 30 kpc. The velocity at an infinite distance 
is either 150 or 220 km s _1 , each being simulated for a 
prograde and a retrograde orbit, thus giving four mod- 
els. The resolution and softening are the same as in the 
unstable disk models 

Our second type of merger model includes several satel- 
lite galaxies taken from a large-scale cosmological simu- 
lation (see Martig et al. 2009a, b). The primary galaxy 
in this case uses 4.7xl0 6 stellar particles, 3.2xl0 6 gas 
particles, and 5.7xl0 6 dark matter particles. The soft- 
ening length for force calculations is 150 pc in this model, 
larger than in the other simulations, but the final thick 
disk scaleheight is still well resolved. The disk is initially 
thin at redshift z = 1.2, and it contains 30% of its mass 
in the form of gas. Over time, the galaxy interacts with 
several dwarf satellites, the most massive of which have 
mass ratios of 11:1, 13:1, 15:1, and 18:1. The final stellar 
mass in the primary is 8xl0 10 M . We analyze its thick 
disk at z = 0. 

4. RESULTS 

4.1. Thick disks from instabilities and giant clumps at 
high redshift 

The initial Gyr in our models produces giant star- 
forming clumps with masses up to a few 10 8 Mq. The 
clumps eventually merge into a bulge (see also BEE07). 
A thick disk forms quickly because of the same gravita- 
tional instabilities that form the clumps. That is, the 




Fig. 3. — Vertical line-of-sight stellar density profiles for the 
systems shown on Figure 2, at r=5 kpc (see Fig. 4 for radial vari- 
ations), normalized to the mid-plane density of the first snapshots 
and plotted in log scale. The initial thin but unstable disk is rapidly 
stirred and thickened by gravitational instabilities and clumps. A 
thick disk is formed within a Gyr, and a faint thin disk forms inside 
of this from left-over gas. Continued external accretion causes a 
more prominent young thin disk to form, and the thick disk shrinks 
a little when this happens. 

instabilities scatter stars and gas to high velocity disper- 
sions, and the clumps that are formed by the instabilities 
continue to scatter stars and increase the dispersion for 
as long as they are present. 

The exponential scale-heights in the stellar disks were 
measured after the clumpy phase, when the clumps have 
either dissolved or coalesced in a bulge. The scale-height 
is ~ 2 to 2.5 kpc. It is constant with radius (Fig. [2]), 
having only a minor flare in the outer regions. That is, 
at a radius of 4 times the disk radial scale-length, the 
scale-height increases to ~3 kpc. 

The slow addition of gas after 1 Gyr causes the 
thick disk to shrink toward the mid-plane, as expected 
(jElmegreen fc Elmegreenl I^OOGa) . In our models, the 
scale-height of the thick disk decreases from ~ 2 or 2.5 
kpc to ~ 1.5 kpc (Fig. 3 and 4). In all models, the thick 
disk remains a separate component and does not become 
a low-density tail to the thin disk (Fig. 3). The thick 
disk also preserves its constant scale-height with radius 
during the accretion (Fig. 4). 

4.2. Thick disks from minor mergers 

The radial profiles of the thick-disk scale-heights in 
our minor merger models are also shown in Fig. 3J 
They have no large thickening in the inner regions (r~0) 
and the outer disk becomes extremely thick, with scale- 
heights of 5 kpc and more. Minor mergers produce 
flared disks instead of uniformly thick disks. This prop- 
erty is not exclusive to our code or initial conditions: 
other published simulations of disk thickening by mi- 
nor mergers show the same thing (I Walker et all ll99l 
IVillalobos fc Helmil l2008t iRead et al.ll2009D . The strong 
flarin g often results in boxy iso photes for edge-on thick 
disks IVillalobos fc Helnfi|[200cl . 

Flaring is inevitable in merger models because the 
dense inner regions have stronger disk forcing compared 
to companion forcing and the perpendicular response is 
small. In contrast, low-density outer regions have lower 
local disk forcing and stronger companion forcing, so the 
disk responds more. Even if an interaction provides the 
same amount of kinematic heating to the whole disk, the 
resultant scale height will vary as the inverse of the sur- 
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Fig. 4. — Radial variations of the thick disk exponential scale- 
height h z are shown as a function of the radius r in units of the 
radial exponential scale-length h r . On the top, models of thick 
disk formation by internal evolution are shown after 1 Gyr (dashed 
lines) and after subsequent thin disk accretion (dotted lines). On 
the bottom, models of thick disk formation by 1:10 mergers are 
shown in four cases (broken lines) and in a cosmological case (solid 
line). The scale height is about constant for models with internal 
disk evolution, but it rapidly increases with radius for models with 
mergers. 

face density. Minor mergers would also add more energy 
and mass to the outer disk than the inner disk because of 
the disruption of the small galaxy. Self-driven disk thick- 
ening by internal clumps does not produce much flaring 
because the perpendicular forcing and the disk gravity 
that resists this forcing are both proportional to the sur- 
face density. 

5. COMPARISON TO OBSERVED THICK DISKS 

Observed thick disks have a relatively constant scale 
height from the central regions to the outskirts of a 
galaxy. Minor merger models do not have this property, 
and they also lack any thick component in the central 
regions. Models where thick disks form by gravitational 
instabilities and clump scattering do have this property, 
however. Therefore, the disky shapes of thick disks are 
explained exclusively by the instability model. 

The scale-heights of our thick disk models are consis- 
tent with the observed values of 1 to 2 kpc for typical 
bright spirals, even after the thick disk shrinks from the 
growth of an underlying thin disk. For Milky Way galaxy 
masses, our models end up with about one-third of the 
initial baryons in a bulge, one-third in the thick stel- 
lar disk, and one third in gas that will contribute to 
the subsequent thin disk. After further thin disk ac- 
cretion, the final masses for the bulge, thick disk, and 
thin disk become 2 - 3 x 10 10 M n , 2 - 3 x 10 10 M Q , and 



8 — 9 x 10 10 Mq, respectively. Thus the thick-to-thin 
disk mass ratio is 30%, roughly typical of modern spiral 
galaxies. 

The thick disks formed in our models are relatively 
massive, 2 — 3 x 10 10 Mq, but consistent with many ob- 
served spiral galaxies. A larger set of lower-resolution 
simulations in BEE07 showed that the fraction of the 
clump mass that ends up in the disk can vary largely, 
so our model can account for a wide range of thick disk 
masses. In this case an anti-correlation between the bulge 
mass fraction and the thick disk mass fraction is pre- 
dicted. That is, clump stars that do not scatter into the 
thick disk will end up in the bulge. YD06 do observe that 
late-type spirals (which have low bulge fractions) have 
higher thick-to-thin disk mass ratios than early-type spi- 
rals (which have large bulge fractions). 

Our models suggest that minor mergers are not the 
main formation mechanism for thick disks. Still, they 
occur for real galaxies over a wide range of redshifts, 
and they probably also occurred when the old thick disks 
formed. Their effect on the thick disk seems to have been 
minor though. One problem minor mergers have is that 
the gas fraction of disk galaxies is high at high redshift 
(z > 1; Daddi et al. 2008), and a high gas fraction can 
prevent minor merg ers from thickening the stellar disk 
(jMoster et al J 120091 ) . A problem at low redshift z < 1 is 
that the thick disk is already in plac e, so minor mergers 
cann ot impact its structure much (|Villalobos fc Helmil 
2008). Also, the mass of a young disk can be larger than 
the mass of a minor companion, and the disk is closer to 
its own stars than the companion is, so the instantaneous 
gravitational force on the disk from its own stars is gener- 
ally much larger than the force from a companion. Thus 
internal processes should dominate thick disk formation 
if young disks are clumpy. This is the essential point of 
the present models: the gravitational potential in young 
disks is extremely clumpy and time-variable because of 
violent instabilities and the presence of spiraling massive 
clumps. 

The stellar disks of spiral galaxies do have some outer 
flaring, but it begins only at large rad i i and it affects 
main ly the thin disk (e.g., lAlardl 120001 : iMomanv et al.l 
120061 for the Galaxy). This situation is consistent with 
a picture in which the thick disk forms mostly by in- 
ternal processes at high redshift, and the primary effect 
of interactions is to warp and flare the outer thin disk 
at lower redshift. Stellar streams within the thick disk 
could also result from dwarf companions (jGilmore et al.l 
120021) . 

An additional consideration is that the young clumpy 
phase of a galaxy should form only one thick disk compo- 
nent that has a single large vertical velocity dispersion. 
If minor mergers form the thick disk, then each merger 
could form a distinction component. The absence of any 
significant vertical gradient of the stellar velocity disper- 
sion (Moni Bidin et al. 2009) suggests there is only one 
thick disk component. 

Thick disk formation by internal forcing is also consis- 
tent with the observed elemental abundances. Our mod- 
els indicate that thick disks and bulges form together 
within a time span of ~1 Gyr. This explains their el- 
evate d a/Fe abundances a nd other abundance similar- 
ities (jMelendez et al.l [2008). If minor mergers were the 
main drivers of thick disk formation, then thick disks 
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would continue to grow below z < 1 (as illustrated by 
our cosmological model). In that case, the age range, 
metallicity range, and formation timescale for the thick 
disk would differ from the observations. Other models 
of thick dis k formation by internal evolutio n, like cluster 
disru ption (lKroupall2002f) or rad ial mixing (jRoskar et all 
2008; ISchonrich fc Binne\l2009 ), would also imply an un- 
acceptable continuity of stellar ages in the thick disk. 

6. CONCLUSION 

There is increasing evidence that many disk galaxies at 
high redshift undergo a clumpy and turbulent phase that 
can drive rapid internal evolution. To understand if this 
is the main mode of spiral galaxy formation, remnants of 
this phase should be idenfitied in the oldest components 
of today's galaxies. The two most obvious of these old 
components are the bulge and thick stellar disk. Previous 
models showed how the clumps formed by gravitational 
instabilities could coalesce into a bulge, but bulges can 
form by galaxy mergers as well. Here we showed how the 
clumps, and the instabilities that formed them, can also 
make the thick disk, while minor mergers fail to explain 
some observed properties of thick disks. 



The standard model in which minor mergers stir an 
existing thin disk and deposit stars far from the plane 
has a problem: it inevitably makes flared thick disks be- 
cause the vertical forcing is independent of the local disk 
column density. Such flaring is contrary to observations. 
Our model produces a uniform disk thickness because 
the forcing is proportional to the column density. This 
property is preserved during subsequent thin disk forma- 
tion by continued slow accretion. In addition, our model 
reproduces the mass fractions of thick and thin disks and 
explains, in general terms, the chemical similarities of the 
Milky Way bulge and thick disk. 

The observed properties of thick disks suggest that to- 
day's spiral galaxies went through a clumpy phase with 
rapid internal evolution, and that clumpy disk galaxies 
at high redshift are progenitors of spiral galaxies. 
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